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We systematically study the localization effect in discrete-time quantum walks on a honeycomb 
network and establish the mathematical framework. We focus on the Grover walk first and rigorously 
derive the limit form of the walker’s state, showing it has a certain probability to be localized at 
the starting position. The relationship between localization and the initial coin state is concisely 
represented by a linear map. We also define and calculate the average probability of localization by 
generating random initial states. Further, coin operators varying with positions are considered and 
the sufficient condition for localization is discussed. We also similarly analyze another four-state 
Grover walk. Theoretical predictions are all in accord with numerical simulation results. Finally, our 
results are compared with previous works to demonstrate the unusual trapping effect of quantum 
walks on a honeycomb network, as well as the advantages of our method. 

PACS numbers: 05.40.Fb, 03.67.Lx 


I. INTRODUCTION 

Quantum walks are generally classified into two cat¬ 
egories: discrete-time and continuous-time quantum 
walks [1-4]. In this paper we focus on the former which 
were first introduced by Aharonov et al. in 1993. A 
discrete-time quantum walk usually involves a quantum 
coin and a walker. The coin controls the walker’s move¬ 
ment. Since a quantum coin could be in any super¬ 
position state of up and down, quantum walks exhibit 
many different features when compared to classical ran¬ 
dom walks (see Refs. [5-10]). For example, the quan¬ 
tum walker spreads at a rate linear in time in the one¬ 
dimensional Hadamard walk [10], much faster than the 
square-root rate in the classical counterpart. 

Another uncommon phenomenon in quantum walks 
is localization. For a classical random walk on a line 
or a square network, as well as the one-dimensional 
Hadamard walk, the probability of finding the walker at 
a specific position approaches zero when the number of 
steps is large enough [3, 11]. However, Tregenna et al. 
numerically discovered that in quantum walks on a 2D 
square lattice, the walker under the control of Grover’s 
operator [12] was highly probable to be found at its initial 
position, which was called the walker’s localization [13]. 
Inui et al. then rigorously analyzed this phenomenon 
and owed localization to degenerate eigenvalues of the 
system’s time-evolution operator [14]. Later localization 
was also found in quantum walks on other kinds of lat¬ 
tices [11] or with other coin operators [15-18]. Additional 
time-dependent phase gates could also cause localization 
in a 2D alternate Hadamard walk according to litera¬ 
ture [19]. 

Some quantum walks have been successfully employed 
to design effective quantum search algorithms [20-24], 
whose target is to transform the system to a specific state 
at a probability of 0(1). Thus the localization effect in 
quantum walks shows potential applications in designing 
fast search algorithms. Some simple quantum walk mod¬ 
els have already been implemented by experiments [25- 


31]. As the structure of graphene, the honeycomb lattice 
has been paid much attention to recently. It may have 
application prospects in quantum controlling and com¬ 
puting based on 2D hexagonal materials like the graphene 
and h-BN [32]. Besides, since the honeycomb lattice is 
more complicated than a line or a common 2D square 
lattice, we could expect more different characteristics in 
quantum walks on a honeycomb network. 

This paper is organized as follows: Sec.H is devoted to 
the definition and mathematical framework of quantum 
walks on a honeycomb network. In Sec. HI, we calculate 
the ket in the limit t —> oo and show the existence of 
localization in the Grover walk, which is in accord with 
numerical simulations. Then in Sec. IV, we investigate 
the dependence of localization on the initial state and 
derive a non-zero average localization probability. Quan¬ 
tum walks governed by other coin operators are consid¬ 
ered and the condition for localization is proposed. Fur¬ 
thermore in Sec. V, an extended model with a four-state 
coin is analyzed. In Sec. VI, we compare our results 
to previous models and show the different properties of 
quantum walks on a honeycomb network. 


II. QUANTUM WALKS ON A HONEYCOMB 
NETWORK 

A. The position space Up and the coin space T-Lc 

The quantum walk of a particle is usually described by 
two Hilbert spaces, the coin space He and the position 
space T-Lp. Different from an infinite square lattice or an 
unrestricted line, a honeycomb network is not a kind of 
Bravais lattice, for here are two types of vertices [33, 34]. 
For example, after we translate the network along the 
vector El in FIG. 1, the new lattice is not identical to 
the original one. Besides, orthogonal Cartesian coordi¬ 
nates, e.g. X-Y in FIG. 1, are inconvenient for describing 
the positions of vertices here. Therefore, we divide the 
network into rhombuses and set up an oblique coordi- 
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FIG. 1. Honeycomb network. Black and white points are two 
types of vertices. The distance between two nearest vertices 
is set to i. 


nate system. A vertex is labeled as {x,y,s), where x 
and y stand for coordinates of the lower left corner of the 
rhombus in which the vertex is contained, and s indicates 
the type of the vertex. We make the convention that a 
black (white) vertex is type-0 (type-1), thus any position 
eigenket is denoted as \x,y, s)p, or |r, s)p in short, where 
f represents (x,y). 

On a honeycomb lattice, a particle is able to move 
in three directions, so the dimension of the coin space 
is three. The eigenstates of the coin are labeled as |1), 
|2), and |3), corresponding to movement along (against) 
El, E 2 and A3 from a type-0 (type-1) vertex. Whether 
the walker will move along or against these vectors is 
determined by the type of the vertex. The system’s total 
state vector is denoted by a ket in the composite Hilbert 
space 'Hc®'Hp. 


B. Time-evolution 


A single step of quantum walks starts with applying 
the coin operator C on the coin space. We first choose 
Grover’s operator G as the coin operator. The three- 
dimensional form of G reads [12] 



-12 2 
2-12 
2 2-1 


Then, the particle is moved by the conditional shift op¬ 


erator S, expressed as [22] 

•S' = XI ® 1)p(gs|, (1) 

j,s.r 


where 


Vi = (0,0), V2 = (1,0), V3 = (0,1). 

Afterward the next step starts. During a single step, 
the position type s is changed by the operator S, so the 
particle will move in opposite directions in the following 
step. Therefore this is a kind of “flip-flop” shift [35]. 

The time-evolution operator U of the total system is 
U = S ■ {C (Slip) where Ip is the identity operator acting 
on the position space. Accordingly the system’s state 
jlkt) after t steps is derived by applying U on the initial 
state for t times: 


IT*) = (£/)* I To). (2) 

The state vector |Tt) can be expanded in terms of all 
the coin and position eigenkets as 

l^i) = ® 

j,r,s 

We introduce the vector 

\Mr,s)) = [V't(f,s),V’?(f,s),V’?(r,s)]'^, (3) 

where the superscript T means transpose. The probabil¬ 
ity of finding the particle at the position (r, s) equals the 
squared norm of s)) if its position is measured. 

Taking into consideration the topology of the honey¬ 
comb network, we can obtain the recursive relations 
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\i^t+i{r,s)) ='^Gj\i;t{r- (-l)®Uj, s © 1)) (4) 

where Gj is a matrix whose jth row is the same as the 
jth row of the coin operator C and other elements are all 
zeros, e.g. 

, r 0 0 o' 

G2 = - 2-12 . 

[0 0 0 

Applying discrete Fourier transformation (DFT, see Ap¬ 
pendix A) to two sides of Eq. (4), we get the recursive 
relations in Fourier space: 

\ipt+i{k,l,0)) = MC\'ij}t{k,l,l)), (5a) 

\^t+i{k,l,l)) = M^C\Uk,lMi (5b) 

where M = Diag{l, and k and I are wave num¬ 

bers corresponding to coordinates x and y. 

Assume the particle starts from position (0,0, 0) and 
denote the initial state of the coin a|l) + /3|2) + qjS) as 
. The initial state of the system is therefore 

|To) = [a,/3,7]'^ ® |0,0,0)p, 


( 6 ) 
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i.e., 

Wo{x,y,s)) = (5s,o<5j/,o'Jx,o[a./3,7]'^.Vx, j/ G Z. 

5a,b is the Kronecker’s delta. After DFT, 

IV'o(ft,Z,s)) = 5 s,o[a,/ 3 , 7 ]'^,VA:,l G (- 7 r, 7 r]. (7) 

According to Eq. (5), 


when t is even. Thus we need to explore the properties 
of the evolution operator tj = MCM^C. 

III. LIMIT OF THE STATE VECTOR 

Since U is unitary, its eigenvalues can be denoted as 
. It is not straightforward to show that = 0 and 
^2 = —^3 = S where 9 satisfies 

cosd = i [4 cos A: + 4cosZ + 4cos(A: — ?) — 3]. (9) 


Eq. (9) essentially agrees with the results in Refs. [36, 37]. 
1,0)) = /, 0)) ( 8 ) The corresponding eigenvectors are 

_I 


I'/'j) 


I 


l^-ik _ 1 

—/+^j) _j_ 2.g*C^j~0 _|_ — g~^^ -j- —g^^j — 






( 10 ) 


where Nj is a normalization factor. 

To determine whether the localization will occur, one 
needs to calculate Poo( 0 , 0 , 0 ), i.e., the probability for the 
particle to be found at its starting point ( 0 , 0 , 0 ) in the 
limit t —>■ oo: 


Poo(0,0,0)= lim(V^t(0,0,0)|V’t(0,0,0)). (II) 

t—^oo 

Obviously, Pt(0, 0,0) = 0,V odd t. As a result, we only 
need to pay attention to even-t cases. 

Substitute the results above into Eq. ( 8 ), 
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( 12 ) 

i=i 

Then the state ket in real space is derived by applying 
inverse DFT to Eq. (12), 


= dMle™*'|^t(A:,l,0)), (13) 


(27r)2 


IV 2 


specially, 


IV’i 


( 0 , 0 , 0 )) = 7 ;^ ^ / dkdle^^^*/^\(l)j){(j)j\Mk,l,0)). 

1^1 dV2 

(14) 

To further simplify the expression of ll/'t(0,0,0)), we 
first introduce a mathematical lemma. 


Lemma 1. Suppose ip{x,y) G L^{[a,b] x [c,dj) and it 
satisfyies that a.e.x G [a, 6 ], ipxiy) = ^{x,y) is a count¬ 
ably piecewise monotonic function of y, and in each 
monotonic interval (c„,(i„), ipxiy) G C‘j^{{cn,dn)) and 
the set {y\dyip{x,y) = 0} is countable. Let f{x,y) = 
g{x, y) + ih{x, y) and g,h € L^{[a, 5] x [c, d]). Then 


lim 

t—^OO 



fd 

/ f{x,y)e^^^^'y'^*dydx = Q. 

J C 


The explanations of some notations as well as the 
detailed proof are given in Appendix B. According 
to Lemma 1, we deduce that the contributions to 
|^/>t(0,0,0)) from items with j = 2,3 in Eq. (14) are 
negligible when t is large enough. This is not hard to 
understand because when t is very large, will os¬ 
cillate much faster than |V’o(A;, ^ 0 )) as fj changes. 

As a result, the integral will become zero when t tends 
to infinity and we only need to consider the term with 
j = 1. Since = 0, 

IV'oo(0,0,0))~/ dMI|</)i)(<^i|^o(A:,/,0)). (15) 
(271-)^ Jv2 

In Eq. (15), only \(fi) is a function of k,l while 
\'tfo{k,l,0)) is constant according to Eq. (7). Conse¬ 
quently Eq. (15) could be interpreted as a linear map 
from the initial state \tpo{k,l,0)) to |^/>oo(0, 0,0)). The 
map is represented by a transformation matrix P 

^= 7 ^ 2 / dkdl\<f,){<f^\. (16) 

The matrix F contains all the information about the 
walker’s behavior at its starting point when t —>■ 00 . The 
existence of localization is directly related to the system’s 
initial state via the matrix P. The final form of F is ob¬ 
tained by exploiting the results in Eq. (10): 



2 -1 -1 

-1 2 -1 

-1 -1 2 


(17) 


Now we take a numerical example. Assuming [a, (3, 7 ]"*" = 
[ 1 , 0 , 0 ]"'", which means the initial state of the coin is 
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“moving along Ei ”, the probability of finding the par¬ 
ticle at the position ( 0 , 0 , 0 ) after a great number of steps 
is |F[1,0,0]'^P = 1/6 accordingly. Apparently, the prob¬ 
ability is not zero and the localization is revealed. 

We also conduct numerical simulations and the results 
are illustrated in FIG. 2. A noticeable peak arises in 
the left of FIG. 2. The right figure manifests the prob¬ 
ability Pt( 0 , 0 , 0 ) oscillates around our theoretical value 
Poo = 1/6 and clearly exhibits a tendency to converge, 
indicating the walker does have a nonzero probability to 
be localized at its starting point. 


IV. CONDITIONS FOR LOCALIZATION 

A. Dependence on the initial state 

We have seen that the limit behavior of the Grover 
walk is represented by Eq. (15). The probability 
Poo ( 0 , 0 , 0 ) is dependent on the coin’s initial state 
[a,/ 3 , 7 ]"'" and it may decrease to zero under some cir¬ 
cumstances. To verify this point, we work out the eigen¬ 
vectors of F: 



The corresponding eigenvalues are 0, 1/2 and 1/2 respec¬ 
tively. So, if the initial coin state is -^[1,1,1]"'", the uni¬ 
form superposition state of the three eigenstates without 
phase differences, named the Grover state [17], localiza¬ 
tion will disappear. Numerical results are illustrated in 
the left of FIG. 3 and Pt clearly converges to zero. It also 
manifests the Grover state is the sole initial state caus¬ 
ing the walker to be delocalized. This result is consistent 
with the conclusion in Ref. [37]. 

On the other hand, the maximal value of Poo (0,0,0) 
is |l/2p = 25%. Any initial coin state which is a linear 
combination of -^[— 1 , 0 , 1 ]"'" and ■^[—I, 2 , — 1 ]"’’ will lead 
to the maximum. For example, when the initial state 
= -^[I, Poo(0,0, 0) will reach 

the maximal probability 25% of localization (FIG. 3). 
This is totally different from the Grover state above. 

B. Average localization probability 

Since the probability Poo(0,0,0) is dependent on the 
initial state of the coin, what is the average probability if 
the initial state is generated randomly? A random pure 
state of the coin can be obtained by applying a random 
unitary operator P 3 G SU(3) on an arbitrary coin state, 
such as 13) (see Sec. IIA). The average probability is 
reasonably defined as 

y 7"oo;a,/3.7dn3 = ^/ 


where dfla is the Haar measure of the group SU(3) and 
fis = f dfla, standing for the volume of SU(3). A unitary 
operator P 3 is usually expressed with the generators Aj 
of SU(3) and the angular parameters (j as [38] 

JJ^ — g*'^3Clgi''‘2C2giA3C3g*'^5C4giA3C5gi''‘2C6g*'^3C7giA8C8 

where 0 Ci, Cs < tt, 0 < C 3 ,C 7 < 27r, 0 ^ C 2 ,C 4 ,C 6 < f 
and 0 ^ (/s ^ \/37r. Under this representation, the Haar 
measure dH 3 of SU(3) and the random initial state PsiS) 
read [39] 
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dH 3 = sin 2(2 cos C 4 sin^ C 4 sin 2 C 6 HdO- 

4 = 1 

gdCi+C 3 ) gQg (^2 sin ^4 

P 3 I 3 ) = —e“dCi-C3) sin(/4 . (19) 

COSC4 

Thus the average probability is obtained by finishing the 
integral in Eq. (18). The result is 



The result manifests that although a Grover state coin 
will delocalize the quantum walker, the expected proba¬ 
bility Poo of localization for a random initial state l^io) is 
still up to 16.7% if we measure the walkers position after 
even-t steps. Thus the localization phenomenon happens 
very frequently in our model and it could be useful in 
designing some quantum algorithms. 

C. Dependence on coin operator 

Apart from the initial coin state, it is obvious that 
the crucial point for the localization to occur is the ex¬ 
istence of highly degenerate eigenvalue 1 of the operator 
U, which is in agreement with literature [11, 14]. This 
characteristic is undoubtedly determined by Grover’s op¬ 
erator G. If there is another coin operator C leading to 
the eigenvalue 1 of U, the other two eigenvalues should 
also be like the form where 6 * is a function of k and 1. 
According to Lemma 1, these two terms in Eq. (14) usu¬ 
ally contribute less to the system’s limit behavior, and 
the term with eigenvalue 1 may give rise to localization. 
Then what are the requirements for such a coin operator? 
We find a sufficient condition. 

Eirst, we extend the model of quantum walks on a 
honeycomb network in Sec.II. In the Grover walk, the 
operator C is independent of position. We now release 
this restriction and let the coin operator be related to the 
position type. The time-evolution operator becomes 

[/ = 5 • (G G ^ |r, 0)(f, 0| + P G ^ |r, l)(r, 1|). (20) 

r r 

Gonsequently, the evolution operator U in Fourier space 
is modihed. 


|PP3|3)|2dL!3, (18) 


U = MDM'^C. 


( 21 ) 
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FIG. 2. (Color online) Numerical simulation results for the Grover walk on a honeycomb network with the initial state 
a = 1,7 = /3 = 0. The left fignre shows the particle’s spatial probability distribution on the honeycomb network after 400 
steps. Points in the right figure are the probability for finding the walker at its starting point as a function of time t and the 
horizontal line is the theoretical resnlt Poo = 1/6. 
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FIG. 3. P-min VS Pmax ■ In the left fignre, localization disappears when the initial state is the Grover state and in the right 
fignre Poo(0,0,0) reaches the maximal value with the initial coin state [1, 


The sufficient condition: li D = and the eigen¬ 
vectors of C are all real vectors (overall phases are omit¬ 
ted), 1 will be one of t/’s eigenvalues and the localization 
is still present. 

The proof is given in Appendix C. The Grover walk 
we have studied above is a special case since C = G = 
= H. Here we also give another example of the coin 
operator 


Specially if we additionally require all the eigenvalues 
of G to be ±1, then we will get the sufficient condition 
for localization when the coin operator is independent of 
positions. 

V. FOUR-STATE QUANTUM WALKS ON A 
HONEYCOMB NETWORK 


H = 


73 


1 1 1 

1 

• 4 -2 

1 


( 22 ) 


H looks much more like the three-dimensional version 
of Hadamard matrix and H satisfies all the conditions 
above. We take C = H and D = with an initial coin 
state [1,0,0]"'" and perform numerical simulations. The 
result is given in FIG. 4, which obviously demonstrates 
that the walker has a nonzero probability to stay at its 
starting point when t > oo. 


A. Definition and transformation matrices 

In this section, we consider another modified model in 
which the dimension of the coin space is enlarged to four. 
The additional fourth component represents not hopping, 
which means the walker is allowed to stay at its current 
position without moving. As a result, the shift operator 
now reads 

3 

S = |4)(4|(g)Ip+^^ |j)(j| (g) |f- {-iyvj,s®l)p{f,s\, 

j—1 s,r 
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and the revised recursive relation is given as 

3 

\ipt+i{r, s)) = C4\il^tir,s))+'^Cj\'tptir-{-l)'''vj,s®l)), 

i=i 


where the coin operator C now should be the four¬ 
dimensional Grover’s operator G [12] 



-11 1 1 ■ 
1-111 
11-11 
. 1 1 1 - 1 . 


(23) 


when the time t is large enough, we only need to pay at¬ 
tention to eigenvalues ±1 and their corresponding eigen¬ 
vectors. The eigenvectors are 

\(j)[) =[e*' - e*^ 1 - e*', -1 -h e*^ 0 , 

e*'=(-l-f e*'),-e*'(-l + e*'=), 0 ]'^ 

corresponding to 1 and 

|.^') = [- 1 , 0 , 0 , 1 ,- 1 , 0 , 0 , 1 ]T 

|,^') = [ 0 ,e-*^ 0 ,-e-*^ 0 ,l, 0 ,-l]T 
|<^;) = [ 0 , 0 ,e-*',-e-*', 0 , 0 ,l,-l]T 


The fourth coin basis state |4) stands for “not moving”. 
Similarly, the recursive relations in Fourier space Eq. (5) 
also change to 

I^t+i(fc, 1, 0)) = C4\Mk, 1, 0)) + M'C\Mk, 1, 1)), 
I^t+i(fc, 1, 1)) = CilMk, 1, 1)) + 1,0)), 

with M' = Diag{l, 0}, which can be further 

expressed as 


corresponding to — 1 . 

Now, if the walker still starts from the position (0, 0,0) 
with an initial coin state [a,/3, 7 ,/x]"'", the initial state of 
the overall system reads 

1 ^ 0 ) = [a,/3,7,M]'^ G | 0 , 0 , 0 )p 
and in Fourier space 

\^o{k,l,s)) = 6s,o[a,(3,'y,fi]'^yk,l G (-7r,7r]. 


\tjjt+iikj,0)) 


Ca M'C' 

'\MkJ,0)) 

JV)t+i(fc,l, 1 ))_ 


M'tCt Ca 

jtptikj, 1 ))_ 


Since U has two highly degenerate eigenvalues ±1, we 
denote 


Therefore, 


'\Mk,l,0)) 

j^tik,l,l))_ 


(ur 


\'ipo{k,l,0)) 

\'^o{k,l,l)) 


(24) 


where U' = 


is an 8 X 8 matrix. So we 


Ci M'C 
M'C"^ Ca 

start to examine the properties of the evolution operator 

ii'. 


The eigenvalues of U' are 1,-1, and where 
01 and 02 are functions of k and 1. According to Lemma 1, 


= dfcdl|<).i)(</)i|, (25) 

(27r)^ Jv^ 

1 f 

^ E )(</., I, (26) 

where vectors are orthogonalized and normalized 
forms of |(()f). Hence, analogous to Eq. (15), for a large t 
we have 


IV’* ( 0 , 0 , 0 ))' 

IV’*(o,o,i)) 


(J-+ -b (-i)‘J-_) 


|'0o(fc,l,O)) 

0 


(27) 
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or 


I^t(0,0,0)) ^ (F+ + (-l)‘F_)|^o(fc,^0)), (28) 

where the elements of 4 x 4 matrices F± are 

{F±)ij = {F±)ij, 1 ^ i, j ^ 4. 

Different from the three-state walk, the extra freedom 
causes the probability Pt(0,0,0) nonzero when t is odd. 
However, Eq. (28) implies the transformation matrices 
for odd and even t are different, due to the negative eigen¬ 
value — 1, i.e. 

i"oo.e(0,0,0)= |Fe|^o(A:,/,0))p (29) 

Poo.o(0,0,0) = |F„|^o(A:,/,0))|2 (30) 

where Fe = F+ + E_ and Fq = E+ — F_. Subscripts ‘e’ 
and ‘o’ stand for ‘even’ and ‘odd’ respectively. 

On account of the complexity of Eq. (26), we only get 
a numerical result of F_: 


F_ 


■ 0.320303 
-0.070303 
-0.070303 
.-0.179697 


-0.070303 

0.320303 

-0.070303 

-0.179697 


-0.070303 

-0.070303 

0.320303 

-0.179697 


-0.17969r 
-0.179697 
-0.179697 ’ 
0.539090 . 


while the accurate form of F|_ is obtained 


F+ 


1 

12 


'2 -1-1 O' 
- 12-10 
- 1-120 
.0 0 0 0 . 


(31) 


As an example, with the initial coin state set to 
^[1,1,0,0]'^, the value of Pt(0,0,0) would converge to 


two values 



Poo,e=[Pe^[l,l,0,0f[^ 

0.222901 

(32) 

Poo.o=[Po^[l,l,0,0]^|2; 

^ 0.092699 

(33) 

according to the derivations above. 

To confirm the re- 


suits, numerical simulations are performed, illustrated in 
Fig. (5). The data points do distribute in two branches 
and get close to our analytical predictions (the two hori¬ 
zontal lines) as the step number t gets large. 


B. Maximal, minimal and average probability of 
localization 

Once we have obtained transformation matrices, the 
properties of this quantum walk are easily analyzed. The 
eigenvalues of Fe are 0.718787, 0.640606, 0.640606, and 
zero. Those of Fg are —0.718787, —0.140606, —0.140606 
and zero, intimating that Pg and Pg have equal maximal 
values, about | ± 0.718787^ « 0.516655, as well as the 
same minimum zero. The eigenvectors of Fg and Fg cor¬ 
responding to zero are both ^[1,1,1,1]"'". So the Grover 
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FIG. 5. (Color online) Simulation results for four-state 
quantum walks with initial state [ 1 , 1 , 0 , (black) and 
|[ 1 , 1 , 1 , 1 ]"^ (blue). 
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FIG. 6 . When the initial coin state is a = /9 = 7 = — 

and ^ the values of Pe and Pg converge to the same 

limit, about 0.516655. 


state is also the sole initial coin state that will delocalize 
the particle in four-state walks, which is shown in Fig. 5. 

What is more interesting is that the eigenvector of Fg 
corresponding to 0.718787 is also the same as that of 
Fg corresponding to —0.718787, whose numerical value 
reads 


If 


max) -[-0.288675, -0.288675, -0.288675,0.866025] ^ 
1 


2 V 3 


[-1,-1,-1,3]' 


This suggests if the coin state is initialized with \(pmax), 
Pg and Pg will reach the maximal value simultaneously. 
The numerical results in Fig. 6 do show that the prob¬ 
abilities of localization for even and odd times converge 
to a single limit, just as we have calculated. However, 
the ket jlft) still oscillates, for Fg\ipmax) is different from 

Fg\^rnax) • 

The maximum of Poo here is around twice as much as 
that for the previous three-state coin model, indicating 
the localization is much stronger. We consider the aver¬ 
age probability again to make a more appropriate com¬ 
parison. In this case, arbitrary unitary operators U 4 € 
















SU(4) and the Haar measure d 04 [38] of SU(4) are in¬ 
volved: 

K = |i^eC/4|4)|2dr!4 « 0.334352, (34) 

p;;=-^ [ |^;C4|4)pd04 « 0.139049. (35) 

“4 J 

Pe is much larger than P = 1/6 in Sec. IV B and Pq 
is only slightly smaller than P. In fact the probability 
Pt{0, 0, 0) is always zero in the previous model when t is 
odd. Thus the localization in four-state quantum walks 
is much stronger than that in the previous model due to 
the additional fourth coin state of not hopping. 


VI. COMPARISONS WITH PREVIOUS WORKS 


Now we compare our results of quantum walks on a 
honeycomb network to other kinds of quantum walks 
studied in the literature. 

The first is work [II] by Inui et al about a particle 
walking on a line under the control of Grover’s opera¬ 
tor. Besides walking left and right, the particle has an 
additional freedom to stay at its current position. So 
the coin space is also three-dimensional. For the situa¬ 
tion considered by Inui et al, we derive the corresponding 
transformation matrix F' 

^ r 1 VG-2 276-5' 

F'= 76-2 76-2 V^-2 . (36) 

[276-5 76-2 1 

Likewise, we calculate the average probability of this walk 
and get the result P' « 0.1684, which is slightly larger 
than 1 /6 for the three-state Grover walk on a honeycomb 
network. However, we notice that the position space here 
is one dimensional but the honeycomb network is a 2D 
lattice. Meanwhile, in our case the walker does not have 
the extra freedom of “stay”. If the walker on a honey¬ 
comb network is also allowed not to hop, just as we cal¬ 
culated in Sec. V, the average probability will get much 
larger. Thus we can say the degree of localization is rela¬ 
tively high in the Grover walk on a honeycomb network. 
The particle can be found at its starting position with a 
good chance. 

On the other hand, we can find out the maximal value 
of ID three-state walk by eigendecomposition 

of F'. The eigenvalues of F' are zero, 76 — 2 and 3 — 76 
and the corresponding eigenvectors are 



Hence the initial coin state ^[1,-2,!]"'’ will delocalize 
the walker, which agrees with the result in Ref. [II]. 
However, the maximum of should be j3 — 76]^ ~ 
0.303 corresponding to the initial coin state -^[1,1,1]^, 
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instead of the value 0.202 claimed in Ref. [II]. Actually, 
one can also get the maximal value 0.303 by submitting 
a = j3 = j = 1/73 into Eq. (8) of Ref. [11]. It is obvious 
that our mathematical method with the transformation 
matrix F is very convenient in studying localization phe¬ 
nomena in quantum walks. 

Apart from the localization probability, another re¬ 
markable difference is that in ID three-state walks, the 
limit state converges to a static wave function, while on a 
honeycomb network, the limit ket oscillates between two 
values. Furthermore, the oscillation in four-state walks 
is not so palpable as that of three-state walks on the 
honeycomb network, since it results from the negative 
eigenvalue —1 of U', the evolution operator in Fourier 
space. 

Kollar et al. in their previous work [40] investigated 
another model where the particle walks on a triangular 
lattice with a three-dimensional coin operator, similar 
to three-state quantum walks in this paper. The au¬ 
thors found when the coin operator is set to Grover’s 
operator, the walker still has a rapidly decaying prob¬ 
ability to be localized at the origin. Analyses suggest 
that coin operators leading to localization will transform 
quantum walks on a 2D plane to quantum walks on a 
quasi-one-dimensional line. In contrast, Grover walks on 
a honeycomb network undoubtedly allow the appearance 
of localization, and we have also found many more coin 
operators leading to walker’s localization in the previ¬ 
ous section. Besides, by comparing FIG. 4 with FIG. 2, 
we can see that walker under the control of FI displays 
different probability distribution from G. There are non¬ 
trivial phenomena present in these models. So we can see 
quantum walks on a honeycomb network do have more 
different features. 

The works [37] by Y. Higuchi et al. and [36] by 
Machida also analyzed the localization in quantum walks 
on a hexagonal lattice. The coin operators are Grover’s 
operator G and G', respectively, in their articles, where 

— 1 — cos e 72 sin e 1 — cos e 
C' = - 72 sin e 2 cos e 72 sine (37) 

[ 1 — cos e 72 sin e — 1 — cos e_ 

and the parameter e G [0,27r). In fact when e is set 
to arcsin(—1/3), C' equals Grover’s operator. Their re¬ 
sults of long-time limit of the localization probability are 
congruent with ours. However, the method of oblique co¬ 
ordinate system used to label the vertices in our article 
is more intuitive and concise. We have conducted more 
detailed analysis such as the average probability and the 
dependence on the coin operator. Moreover, we can see 
that the eigenvectors of C" are 

— 1 cot e -I- CSC e I — cos e 

0 , —72 and 72 sine , (38) 

I cot e -I- CSC e I — cos e 

and G' = G'\ So C satisfies the sufficient condition in 
Sec. IV G. Thus I is one of the the time-evolution oper¬ 
ator’s eigenvalue and localization is probable to occur in 
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FIG. 7. (Color online) Relationship between average position 
f and time t. Black and bine dots stand for Fig. 2 and Fig. 4, 
respectively. 

Machida’s model according to the more general result in 
Sec. IV C. 

As a by-product, we numerically estimate the spread¬ 
ing speed of the walker mentioned in their papers. Rea¬ 
sonably, we use the average radius f defined as 

^ = X! y, s) (39) 

x,y,s 

from the starting point to describe how far the walker 
has reached, where 

rx,v,s = y^(x-b|)2-b(y-b^-^)2 

is the distance between position {x,y,s) and (0,0,0). 
The results for the circumstances of Fig. 2 and Fig. 4 are 
plotted in Fig. 7, where the data points are distributed 
in two lines very well, i.e., the walker is linearly spread¬ 
ing on the network and the radius of the ring in Fig. 2 
and in Fig. 4 are proportional to time, which agrees with 
Refs. [36] and [37]. Additionally, the difference between 
the slopes of the two lines indicates the spreading speed 
is related to the coin operator. 

VII. SUMMARY 

In this paper, we have established a precise mathemat¬ 
ical description of quantum walks on a honeycomb net¬ 
work and conducted a comprehensive and rigorous study 
on the localization effect. We first set the coin operator to 
Grover’s operator. With discrete Fourier transformation, 
we analytically obtained the system’s state vector. The 
localization probability in the limit of large number of 
steps is connected to the initial coin state through a con¬ 
cise transformation matrix F, which completely contains 
all the information about the walker’s limit behavior. We 
then analytically showed localization and delocalization. 
The results were all supported by numerical simulations. 
The average probability of localization’s occurrence is 1 /6 


if the system’s initial state is randomly generated. Other 
coin operators were also discussed and we derived a suf¬ 
ficient condition for localization. Further, we extended 
the coin space to four and conducted parallel calcula¬ 
tions, finding an oscillating limit state and deriving a 
larger probability. Some other models in literature were 
reviewed. Basing on the average probabilities, we con¬ 
cluded that the trapping effect of a honeycomb network 
is relatively strong. Our mathematical method with the 
transformation matrix is very convenient to study local¬ 
ization in quantum walks. We numerically showed the 
walker’s linear spreading property as well. 
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Appendix A: discrete Eourier transformation 

Suppose X(rn) is a function of m, an n-dimensional 
vector whose components mj G Z. Its discrete Fourier 
transformation is defined as 

X{k) = (Al) 

rrij GZ 

where k is an n-dimensional wave vector with components 
kj G (—7r,7r]. The inverse transformation is: 



The integration region U„ in Eq. (A2) is kj G (—tt, tt], rep¬ 
resenting an n-dimensional hypercube whose sides equal 
27r in the Fourier space. 

Appendix B: Proof for Lemma 1 

Mathematical notations are as follows: 

• L^([a, 6]): the set of Lebesgue integrable real func¬ 
tions defined on the interval [a, b], 

• Cg{A): the set of bounded real functions defined on 
set A with continuous first and second derivatives. 

• S[ll([a, 6]): the set of Lebesgue measurable real func¬ 
tions defined on the interval [a, 5]. 

• P{x), a.e. X G A: the proposition P(x) does not 
hold only for x G Ai C A and the Lebesgue mea¬ 
sure of Al is zero. 
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Lemma 2 (Riemann-Lebesgue lemma). If f{x) G 
L^{[a,b]), then 

fb 

f (x) cos (tx) dx = lim / f(x)sm(tx)dx = 0. 

J a 



Lemma 3. If a countably piecewise monotonic func¬ 
tion ip{x) definied on [a, b] satisfies that in each mono¬ 
tonic interval {an, bn), ^p{x) G C'^((an,6„)) and the set 
{x I y}'{x) = 0} is countable, then 


lim 

t—¥CX) 


/ cos {ipix)t) dx = lim 

/ '.-r \ , 

a 


sin {Lp{x)t) dx = 0. 

(Bl) 


Proof. We prove it in the case of cosine. The domain of 
_I 


ip{x) can be divided into countable intervals, i.e., 

[a, b) = bb) and [a„, 6„)n[am, bm) = 0, for n ^ m, 

n 

and in each interval the function is monotonic. In this 
case, each piece is bounded and differentiable to the sec¬ 
ond order. 

Consider one of those intervals [a„,&„). First, based 
on the theorem of existence of inverse function, there is 
an inverse function of ip, called (p~^ ■ Besides, near the 
end points a„ and bn, we have 

lim / |cos ((/3(a:)f)| dx = lim / |cos ((/9(x)t) | dx = 0, 


where a'n = an -h e and b'n = bn — e. Because the set 
{x I (p'{x) = 0} is countable, p'ia'n) or ip'{bn) is not zero 
for almost all e. Now in the interval [a'n, b'n), by making 
a substitution u = p{x) and integrating by parts, we get 


cos {ip{x)t) dx = 


cos (ut) , 1 sin (ut) 

du = 


v«) P'Vp~^W\ ip'[v~^{u)\ 


I (sin{ip{b'n)t) sin{p{a'n)t) 

t 


P'iK) 




viK) I MK) ip"[ip-^{u)] , 

+ -- sin{ut) ^ ^ [ du 

p''[p~Hu)] , 

/ smiut) -^- ^du. 

’{p'[p-^{u)]f 


(B2) 


r 


Since |sin(ut)| < 1 and p'{a'n) ^ d,(p'{b'n) ^ 0, the first 
term tends zero when t —>■ oo. For the second term, 
making inverse substitution again, 




h«) W[v Hw)]} 


P''{x) 

W{x)f 

1 1 


dx 


(B3) 


v'iO v'iKY 


According to Lemma 2, the second term tends to zero 
when f —>■ oo. 

Now Let us see the integral in interval [a„, 6„) 

r^n r^'n r^'n 

/ cos {(p{x)t) dx = ( / -\- -\- ) cos {ip{x)t) dx. 

J an J an ^ a'n ^'n 

After letting e —>■ 0 and t —>■ oo, 

/•bn 


lim 

t—¥00 


cos {ip{x)t) dx = 0. 


Finally, since lim |cos (i^(x)t)| dx < b—a is uniformly 
convergent, the limit and the sum is commutable. So, 

pb _ pbn 

lim / cos ((p(x)t) dx = lim ^ / cos ((p(x)t) dx 

J a i—^oo ^ J an 

fbn 

= lim / cos {ip{x)t) dx = 0. 

^f —>-00 / 


□ 

Lemma 4. For a countably piecewise monotonic func¬ 
tion p{x) satisfying that in each monotonic interval 
{an,bn), ip{x) G Cg{{an,bn)) and the .set {x|v3'(x) = 0} 
is countable, if f{x) = g{x)-\-ih{x) and g,h G L"'^{[a,b]), 
then 

nb 

lim / /(x)e*‘^(")‘dx = 0. (B4) 

Proof. 

/(x)e*‘^^®^* =g{x) cos {(p{x)t) — h{x) sin {(p{x)t) 

-\-i[g{x) sin {(f{x)t) -\- ih{x) cos (v3(x)t)]. 

(B5) 

For the first term in Eq. (B5), since g G L^{[a,b]), Ve > 
0, 3 a staircase function gi(x) satisfying 

/ - ffi(x)| dx < e. 

J a 

) pb 

g{x) cos {ip{x)t) dx — / (/i(x) cos ((/?(x)f) dx 

J a 

< / \9{x) - 5 i(x)|dx < e. 

J a 


Hence 
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Because gi is a staircase function, according to Lemma 3, 

nb 

lim / ( 7 i(x) cos ((/9(a;)t) dx = 0. 

t^OO I 
J a 

Therefore, 


where G(x) is integrable. According to Lemma 5, we get 

= lim [ I{x,t) dx = 0. 


pb pd 

lim / / /(x,y)e*‘^(^'2')*dydx 

Jc 


□ 


fb 

lim / g{x) cos {(f{x)t) dx = 0. 

i—^oo / 

a 


Appendix C: Proof for the sufficient condition 


The conclusions are the same for the other three terms 
in Eq. (B5). □ 


Lemma 5 (Lebesgue’s dominated convergence theorem). 
Suppose {fnix)} C 91t([a, 6]) and lim fn{x) = fix), a.e. 

n—^oo 

X € [a, 6]. If fix) is dominated by an integrable function 
gix) in the sense that |/(x)| < gix), then 

pb pb 

lim / fnix)dx= / fix)dx. 
n—^oo I I 

^ a J a 


Further, the index can be continuous, namely, 

j'b pb 

t—¥00 


pb pb 

lim / fix,t)dx= / fix)dx. 
da da 


Denoting the normalized eigenvectors of the matrix C 
as \‘Pm) and corresponding eigenvalues as (m = 

1,2,3), we can construct a unitary matrix P by ar¬ 
ranging three eigenvectors \'Pm) together, namely, P = 
IA 2 ), Ias)]- Then 


Pt/pt = (PMPt)(PGl'pt)(PMtpt)(PC'Pl') 


ip 

'1 0 

0 e-* 

1 

o o 

Pt) 

0 0 

0 0 


o 

o 

e"*' 


0 0 




0 

o' 



0 

0 ■ 

iP 

0 

gik 

0 

pt) 

0 


0 


0 

0 

e*' 


0 

0 



We denote the unitary matrix 


Now we prove Lemma 1: 
Proof. Set 


Iix,t)= r fix,y)P‘^^-’y^*dy. 
J C 

According to Lemma 4, one can find that 
lim /(x,t) = 0, a.e.x G [a, 5]. 


Besides, 


|7(x,t)| = 




< [ \fix,y)\dy = Gix), 

J c _ 


'l 

0 

0 ■ 


-g-ia;i 

0 

0 

0 

g-*k 

0 

pt) 

0 

g-i<V2 

0 

0 

0 



0 

0 

g-i<V3 


as J. Given that eigenvectors \iprn) are all real vectors, 
the matrix P would be a real matrix. Then, PUP^ = 
J J*. The determinant of J J* — I is 

det(JJ* - I) = det(J(J* - J^')) = det(J)det(J* - J'^) 

Since the matrix J* — is an anti-symmetrical three- 
dimensional matrix, its determinant is zero. Therefore, 
we have proved that 1 is one of the eigenvalues of the 
matrix PUP\ as well as the matrix U, because unitary 
transformations will not change the eigenvalues of an op¬ 
erator. 
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